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Abstract 

A general transport equation for the center of mass motion is constructed to study transports of 
electronic system under uniform magnetic field and microwave radiation. The equation is applied 
to study 2D electron system in the limit of weak disorder where negative resistance instability 
is observed when the radiation field is strong enough. A solution of the transport equation with 
spontaneous AC current is proposed to explain the experimentally observed Radiation-Induced 
Zero Resistance State. 



The discovery of the zero-resistance states (ZRS) in two dimensional electron gas under 
uniform magnetic fielc^^ and microwave radiation has triggered a lots of theoretica Pi^i^i^i'^'i^i'^^ 
and experimental activities^iiS to understand the origin of this nontrivial state. Most of 
the theoretical work suggested that the origin of the ZRS is closely related to a negative- 
resistance instability that occurs in the system due to the combined effect of quantized 
Landau levels and photon-assisted scatterin g^i^i^i^i-^^ . It was proposed that the ZRS can be 
explained if the current-dependent resistance of the system which becomes negative at small 
current (for strong enough microwave radiation) becomes positive again when the current 
j becomes large enougb^*^^. The above physics was put together phenomenologically into 
an equation 

E, = ph[jxz] + R[\j\]j (1) 

where R[j] is a phenomenological current-dependent resistance which is negative at j = 0, 
increases as a function of j and passes through zero at \ j\ = jo^, is the applied DC electric 
field and pu is the ordinary Hall resistivity. It was shown that equation (0) admits time- 
independent, stripe-like spatially inhomogeneous solution which leads to zero differential 
resistance for net DC current less than a threshold value^. An obvious theoretical question 
is whether Eq. with the required property of R[j] can be derived microscopically. This 
is the subject of this paper. 

Starting from first principle we shall derive in the following a transport equation for the 
center of mass velocity v = j/ne that treats the effect of radiation to all order with the 
only expansion parameter in the problem being the strength of disorder. We note that a 
transport equation can also be derived from a Quantum Boltzmann Equation approach^i. 
However it is difficult to obtain clear, analytical result in this approach because of the 
intrinsic complexity of the Boltzmann equation formulation itself. We shall see that the 
equation of motion for the center of mass offers a much simpler alternative. 

Our approach to the transport equation begins from the known observation that there 
exists an exact, one-to-one mapping between the solution of the Schrodinger equation of 
a (charged) many-particle system in the absence of the microwave radiation and in the 
presence of the radiation for a class of Hamiltonian^'iiSiiLi^^ 

H{t) = Ho-eJ2ri.E{t) (2) 
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where 

Ho{n) = E ^{P^ - lm)f + \ E Vi^^ - + ^ E r^■K.n (3) 

i ij^j i 

where pi = —ihWi is the canonical momentum for the z*'' particle in the system. A{r) = 
— |(r'x B) is the vector potential corresponding to a uniform, time-dependent magnetic field, 
E{t) is a time- dependent uniform electric field and V{r) is the interaction potential between 
particles. The last term represents an external harmonic potential acting on the particles. 

The physics of the exact mapping can be seen by performing a coordinate transformation 
to the center of mass frame of the many-particle system. In the non-relativistic limit, the 
wavefunctions in the laboratory and CM frames are related by 

^za6(r.; t) = iJcMin - R{ty, t)e'^(^W), 

where R{t) = < > is the center of mass coordinate in the laboratory frame, 

= number of particles and 6{R{t)) is an overall phase that depends on R{t) only. The 
corresponding Hamiltonian in the CM frame is 

H'{t) = E - mhit) - -A'aMm' + i E ^(^1 - (4) 

+ 1 E(^1 + myK-if". + R{t)) - e E ^..E'{t) 

where P = f-R{t),p'- = -ihW and A^f) = -\r'^B'. B' = BandE'{t) = E{t) + \R{t)^B 
are the Galilean transformed magnetic and electric fields in the CM frame valid in the non- 
relativistic limit. The Schrodinger equation in the CM frame, ihipcM = H'(t)ipcM can 
be simplified by a performing a gauge transformation "tpcM = <Pcmg^^^^\ where hA{t) = 
Ei [mR{t).r' -\jUt'R{t').KR{t')^. With Eq. (jH), we obtain 

in^^ct^cMif'-.t) = (H^if^^ - E^1-«w) <PcM{^;t), (5) 



where a{t) = e yE{t) + ^{R{t) x B)j — K.R{t) —mR{t) = 0. a(t) vanishes because it is the 
equation of motion for R{t) derived directly from the corresponding Heisenberg equation 
of motion. We thus arrive at the conclusion that for the class of Hamiltonian 0, there 
exist a one-to-one mapping between the solution of the Schrodinger equation in the presence 
of the microwave radiation %p{ri\t), and in the absence of the radiation (pcMi'f^i'jt), where 
i^i^ut) = 4'CM{^i — R{t);t)e^^^*\ Physically, for the particular form of Hamiltonian we 



considered, the wavefunction of the system follows the center of mass motion rigidity in the 
presence of the radiation field. 

The above result suggests that for more general Hamiltonians of form Hdt) = H{t) + U, 
a new perturbation scheme where the microwave radiation is treated exactly to all order can 
be set up by treating U as perturbation. The perturbation scheme can be set up most easily 
in the center of mass frame. We shall consider static impurity potential U = XJ2iu{ri) in 
the following. Notice that a static potential become time-dependent in the CM frame and 
should be treated by time-dependent perturbation theory. 

To derive the transport equation we start from the exact Heisenberg equation of motion 
for the center of mass coordinate, ihR{t) = [R(t), Hcit)] with K = 0*. We obtain after 
some simple algebra 

mR{t) = e (^E{t) + ^(Rit) x ^)) - ^ / ^'^^ (VM(f)) n(f, t) (6) 

where n{r,t) = < 5{T — ri) > is the time-dependent average electron density. In the CM 
frame where n{f,t) = ncM{r — R{t),t), we obtain for small A from linear-response theory 

ncM(r", t) = n^lin + \ [ rfVrft"x(r-^ - r"" ; t - t'Mf^' + R{t')), 



where xi'f^ ]t ~t') is the (equilibrium) retarded density- density response function in CM 
frame derived with the Hamiltonian Ho- Going back to the laboratory frame and performing 
the disorder-average < u{'r>= and < u{r)u{f') >= |Mp5(r — ■P'), we obtain to second order 
in A, an impurity-averaged equation of motion for R{t) in laboratory framed, 

mR{t) = e (^E{t) + ^{R{t) x B)^ + aV^(^) J dt'x{R{t) - R{t'); t - t') (7) 

where n = N/V and a = X'^\u\'^/n. The equation is manifestly gauge invariant and suggests 
that to second order in the impurity potential, the effects of particle statistics and interaction 
are reflected only in the density-density response function. In the following we shall apply 
this equation to study the ZRS in 2D electron systems. 

In this case we consider E{t) = Eo sin{ujt) + E^, where E^ is a small DC electric field. To 
simplify the equation we divide the center of mass motion into "fast" and "slow" parts, 

R(t) = Rfit) + Rs(t) ~ Ao cos(cjt) + Bo sin(cjt) + Rsit). (8) 

Rf{t) describes the center of mass motion induced by the radiation field whereas Rs{T) 
describes induced motion under the DC field E^. The two kinds of motions are coupled by 



the impurity scattering term which is a non-hnear function of R{t). We shall assume that the 
coupling between Rs(t) and Rf(t) does not modify qualitatively the behavior of Rf(t) and 
the main effect of coupling is to produce an effective equation of motion for Rs{t). Notice 
that we keep only the first harmonic terms in Rf{t) in Eq. (jHJ. This is valid in the small 



Eo limit when the size of the orbit -Rc ~ V^o + -^o much less than the magnetic length 
I = (hc/eBy^"^ and is justifiable under the experimental conditioni*^ where the magnetic 
field is weak and the magnetic length is very long. 

To treat the impurity scattering term we write V^^-^^ /Ig^^ — R(t');t — t') = 

— J qd'^q J ^/mx(g, cu') dt' ^^^■{R{'t)-^R{'t'))-'^' i't-'t')] ^ where cu) is the Fourier transform 
of t)- To evaluate the integral over t' we make the local approximation Rs{t) — Rs{t') ~ 
v{t){t - t') + \^{t){t - and 

/* (H'(^mm-R{t'))-^'{t-t')] ^ /■* ^^/g»[g.(i?^.(t)-i?^(t'))+(q.^7(t)-^'){t-t')] 1 1 I iq.'vjt) ^^ ~ ^'^"^ ^ 

J —oo J—oo \ 2, J 

valid for v{t + T') — v{t) << v{t), where T' ~ 2tt/ujb is the characteristic time-scale for the 
density response. The Rf(t) — Rf{t') term can be expanded in a series of Ao,Bo using the 
identity e'^'^™^ = J2m Jm{x)e^"^^. We obtain after some algebra, 

V«^(,)/_* dt'xim-Rit');t-t') = ^ e^("')(-*+^)|gt/V„(^(g-))J„,(^(?1) (9) 

^ in,m' 



I .q.vit) 92 > 



where z{q) = y {q-AoY + {q.BoY and tan6' = q.Ao/ q.B^. We shall consider zi^q)"^ = R^q'^/2 
in the following corresponding to circularly polarized light^. In this case the response of the 
system is isotropic and the computation is much simplified. The effective force on Rs{t) is 
derived from the m = m' terms in (p)). To order {Rc/lY "we keep only m{m') = 0, 1 terms. 
Putting it back into the equation of motion (jZj), we obtain 

(m + n[a;, \v{t)\]) ij{t) = eEa + ^{v{t) x B) - R[uj, \v{t)\]v{t), (10) 

where ut is the frequency of the microwave radiation, and 

R\uJ, v\ ^ j q.vd'^q 
t>^ J 



ziq)'^ z{q\^ 
1 —)Imx{q, q.v) + -^Imx{q, uj + q.v) 



(11) 



z(q)'^ / ^{q}'^ I 

)Rex{q,uj') 

u)'=q.v 



u)' =u)+q.v 
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It is obvious that R[uj,v\ represents dissipative response of the electron gas to external 
perturbations whereas n[a;,'y] is an effective mass correction on the center of mass motion 
coming from the corresponding reactive response. Negative contributions to the resistance 
shows up in the second term of R for density-density response functions with resonant 
structure, u) ~ 

En SiM^^, where cj„ > are resonant energy levels and 5'„(g) is a 
positive definite function. In this case, R ~ En {Rin^n + R2n{^n ~ ^)), where Rin and i?2n 
are positive definite numbers. We see that negative contributions to R2 exist for uJn < oj. The 
physical origin of the negative resistance has been discussed in several earlier work a^i'^i^i^i-^^ 
and we shall not repeat them here. Correspondingly, the effective mass contribution from 
level n is positive (negative) when uj > {<)ujn- We note that the effective mass correction 
originates from the impurity scattering term is of order n[co', f] ~ m x (l/cu^r) << m in the 
weak-disorder limit, where ujb = eM/mc is the cyclotron frequency and r ~ (fi'2d'^^|'"P)~^ is 
the elastic lifetime. Eq. pOj) differs from the phenomenological equation mainly in the 
presence of the inertial term mv which allows the system to admit time-dependent solutions 
in the present case. 

To see whether R[uj,v] has the expected behavior we consider den- 
sity response function of non-interacting gases where^ u) = 

En.,nK!/ni!) (gV2)"^-^ e-(QV2) (^^^^"^(QVa))' ( "X'-^y:!^."^ ' ^^^^^ 
rii = max(?T,, m), n2 = min(n, m), L''^{x) is the associated Laguerre polynomials, riple) 
is the Fermi distribution function and = To incorporate inelastic lifetime 

effects we also introduce a phenomenological broadening T(T) to the Landau levels, 
i.e. 6{e — uub) (^)^^r/((e — tujJbY + T^). R[uj,v\ is evaluated numerically with 
these approximations. In figure one we present numerical results for the normalized 
resistance R{ijJ,v) / R{Q,Q), as a function of normalized velocity vn = v/{Iujb) for 
uo/ub = 0,0.85,1.0,1.1,1.2, taking {Rc/lf = 0.1, Ep ~ IOub, T = 2ujb, T = 0.5 and 
keeping 20 levels in the sum. We observe that R{uj, v ^ 0) becomes negative when u ^ ujb, 
in agreement with previous results^i. For uj > ujb R{uj,v)^s increases and cross zero at 
around Von ~ 0.05. The effect of microwave radiation decreases rapidly for vn ~ 0.2. These 
qualitative behaviors of R{uj, v) are in agreement with expectation and are not modified by 
changing T or F. 

Eq. (fTUj) allows time-dependent solutions. In the absence of the DC field, a simple, 
spatially homogeneous solution which allows the system to stabilize itself around the point 
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Resistance versus velocity 

for different microwave frequencies 
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1: 



Normalized resistance as a function of velocity for oo/ujb 



0(i), 0.85(m), l.O(m), l.l(w), l.2{v). 



{Iujb)VoN is 



with Ro = Vo/uJs where Us = eB /{m + Ii{uj, Vo))c ~ ujb- The solution represents a collective 
circular motion of the whole fluid moving with speed Vo- Notice that a time- independent, 
spatially inhomogeneous solution corresponding to a pattern of alternating current stripes^ 
still exist. However this solution is energetically less favorable because it requires a higher 
energy to create the charge inhomogeneity needed to maintain the stripes of currents. 

To understand the ZRS we have to consider the boundaries of the droplet of electron 
fluid. For sharp boundaries the boundary condition j± = has to be imposed where j± is 
the component of current perpendicular to the boundary. As a result an edge region with a 
time-independent current j// ~ jo must form. The size of this region is determined by the 
microscopic charge dynamics^ which is still undetermined. Nevertheless according to Eq. 
(fTUj) an electric field perpendicular to the boundary with magnitude ^ = Bvo/ c has to 
be present in this region to maintain the steady current flow. A similar edge region also 
exists at the opposite edge with a current running in the opposite direction, rather similar 
to edge states in Quantum Hall Effect. 

A state with a small net current flow can be created with minimal disturbance to the 
system by shrinking the size of one edge region and enlarging the other. In this case, the 
net voltage drop across the sample is given by 



y ~ 
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corresponding to a resistance matrix with p^x = and pxy = pn, i-e. the ZRS. 

Some comments about the vahdity of our theory is in order. We note that in deriving 
Eq. (jiup we made the local approximation which assumes slowly varying v{t) whereas the 
spontaneous current state we propose oscillate with a frequency ~ ujb- The oscillatory 
solution is allowed because of the presence of the inertia term mv{t). Our analysis shows 
that the correction to the inertia term is small (~ m/{ujBT)) in the limit of weak disorder 
and the local approximation mainly affects R[ljJ, v\. Therefore our general description of the 
ZRS should remain valid as long as the qualitative property of R[uj,v\ is correct. Another 
simplification we employed in our analysis is the assumption of circularly polarized light. The 
response of the system which is isotropic in this approximation would become anisotropic 
when this assumption is relaxed^^. 

Lastly we made a comment on the macroscopic nature of the spontaneous current state 
we proposed. We note that in general a spontaneous current state with \v{t)\ = Vo is 
characterized by a position and time dependent (2D) unit vector field n{x^t) representing 
the direction of the current. The order parameter field n{x, t) has the same symmetry as the 
ordinary 2D x — y model, or superfluids. The main difference between the ZRS state and 
superfluids is that the rigidity of the order parameter is protected by repulsive interaction 
in the case of superfluids, whereas it is protected by the principle of least dissipation in the 
ZRS. The similarity between the two systems suggests that the two systems may share some 
common macroscopic features. For example, vortex-like solitonic excitations may exist in 
the ZRS and may lead to the residue resistance ~ Roexp{—(To/T)) observed in the ZRS 
stated. The existence and nature of solitonic excitations depends on the detailed current 
dynamics of the ZRS state and will be investigated in a coming paper. 
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